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The numerical experiments of turbulence conducted by Gotoh et al. 
are analyzed precisely with the help of the formulae for the scaling ex- 
ponents of velocity structure function and for the probability density 
function of velocity fluctuations. These formulae are derived by the 
present authors with the multifractal aspect based on the statistics that 
are constructed on the generalized measures of entropy, i.e., the exten- 
sive Renyi's or the non-extensive Tsallis' entropy. It is shown, explicitly, 
that there exist two scaling regions, i.e., the upper scaling region with 
larger separations which may correspond to the scaling range observed 
by Gotoh et al., and the lower scaling region with smaller separations 
which is a new scaling region extracted first by the present systematic 
analysis. These scaling regions are divided by a definite length approx- 
imately of the order of the Taylor microscale, which may correspond 
to the crossover length proposed by Gotoh et al. as the low end of the 
scaling range (i.e., the upper scaling region). 

47.27.-i, 47.53.-fn, 47.52.+j, 05.90.+m 

After the discovery of the Kolmogorov spectrum [Q, there have been plenty of 
trials including ||^,^ in order to understand the intermittent evolution of fluid in fully 
developed turbulence. Among them, the line based on a kind of ensemble theoretical 
approaches, initiated by the log- normal model consists of the /3-model (a uni- 

fractal dimensional analysis) 0, the p-model (a multifractal model) the 3D 

binomial Cantor set model and so on. On this line, an investigation of turbulence 
based on the generalized entropy, i.e., Renyi's or Tsallis' IT^JT^ , was started by 
the present authors p^-pO|. After a rather preliminary investigation of the p-model 
[l^ , it has been developed further to derive the analytical expression for the scaling 
exponents of velocity structure function |15| - |l8|], and to determine the probability 
density function (PDF) of velocity fluctuations |18|-|2C|] by a self-consistent statistical 
mechanical approach. 

With the help of the analytical formulae derived in [p"5|-p0[, we will analyze, in 
this paper, the PDF's of velocity fluctuations observed in the beautiful DNS (i.e., the 
direct numerical simulation) conducted by Gotoh et al. [|T]. We will deal with the 
data at the Taylor microscale Reynolds number R\ — 381, since at this Reynolds 
number Gotoh et al. observed the PDF with accuracy up to order of 10^^ ^ 10^^°, in 
contrast with any previous experiments, real or numerical. We showed in [l9|] that our 
formulae can explain quite well the PDF's observed in the real experiment by Lewis 
and Swinney [|2| for turbulent Couette- Taylor flow at R\ = 270 (Re = 5.4 x 10^) 
produced in a concentric cylinder system in which, however, the PDF's were measured 
only with accuracy of order of 10~^. Note that the success of the present theory in 
the analysis of this turbulent Couette- Taylor flow may indicate the robustness of 
singularities associated with velocity gradient even for the case of no inertial range. 
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We have already made clear in that the present theory can also explain quite well 
the PDF's of longitudinal velocity fluctuations reported by Gotoh et al. ||2l[, and have 
revealed the superiority of our PDF to the one derived in [ p3| when the accuracy is 
raised up to the order of 10^^. Assured by the success of these rather preliminary 
tests, we will apply our theory for further precise analyses of the data obtained in 
pljl including the PDF's of transverse velocity fluctuations in addition to those of 
longitudinal fluctuations. 

The basic equation describing fully developed turbulence is the Navier-Stokes equa- 
tion du/dt + (u • V)u — —V {pI p) + v'S/^u of an incompressible fluid, where p, p and 
V represent, respectively, the mass density, the pressure and the kinematic viscosity. 
Our main interest is the correlation of measured time series of the streamwise veloc- 
ity component u, say x-component, of the fluid velocity fleld u in the turbulent flow 
produced by a grid with size £o putting in a laminar flow. Within Taylor's frozen flow 
hypothesis, the quantity of our interest reduces to 5u{r) = \u{x + r) — u(x)\ repre- 
senting the spatial difference of the component u. We assume that in the downstream 
of the grid there appears a cascade of eddies with different sizes = ^n^o where 
5n = ((5 > 1, n = 0, 1, 2, ■ • •). At each step of the cascade, say at the nth step, 
eddies break up into 5 pieces producing an energy cascade with the energy-transfer 
rate e„ that represents the rate of transfer of energy per unit mass from eddies of size 
£„ to those of size In+i- 

The Reynolds number Re of the system is given by Re — 8uqIqIv — {Iq/vi)'^/'^ 
where rj = (z/^/e)^''^ is the Kolmogorov scale ^ with the energy input rate e to the 
largest eddies of size £oj i-e., cq = e. Here, we introduced the notation (5w„ = Su{in) 
representing the velocity difference across a distance r ^ £„. Then, our main interest 
in the following reduces to the fluctuation of velocity difference Sun corresponding to 
the size of the nth eddies in the cascade. Note that the dependence of the number of 
steps n on r/rj, within the analysis where intermittency is not taken into account [Q, 
is given by 

n= -log5r/?7-f (3/4)log5Re. (1) 

For homogeneous and isotropic turbulence, there is a relation between the Taylor 
microscale Reynolds number Rx and the Reynolds number Re |2^ : i?| = ARe where 
A is a real number of the order of 0.01 ^ 10 depending on the experimental setup. 

For high Reynolds number Re ^ 1 or for the situation where effects of the kine- 
matic viscosity can be neglected compared with those of the turbulent viscosity, the 
Navier-Stokes equation for incompressible fluid is invariant under the scale transfor- 
mation (2|J|1: f ^ X°'/^u, \r,t^ \^'°'/^t and (p/p) A^"/^ (p/p). Here, the 
exponent a is an arbitrary real quantity which specifies the degrees of singularity in 
the velocity gradient \du{x) / dx\ = lim^^^o \u{x + £„) — u{x)\/£n = limf„_,o 5un/£n- 
This can be seen with the relation Sun/Suo = {(n/£o)"^^ , which leads to the singu- 
larity in the velocity gradient for a < 3, since bunjGx"^ — const.. We also have 
the relation en/e = (^n/-^o)"^^ ■ 

Let us determine the probability p(")(Q;)(ia to find at a point in physical space 
an eddy of size £„ which has a value of the degree of singularities in the range 
a a + da. Assuming that each step in the cascade is statistically independent, 
i.e., P'^'^\a) = [P*^^^(q;) " our task reduces to determine P^^\a). In order to pro- 
ceed, we have assumed [isl-^ that the underlying statistics describing the intermit- 
tent evolution of fully-developed turbulence is the one given by the Renyi entropy 
ll| S^[P'^^\a)] = (l-g)"Mn/ daP(i)(a)«, or by the Tsallis entropy ill,|ll,|3 

S'^[P^^\a)] — {\ — q) ^ (J daP*^^) (a)"^ — l) . Taking an extremum of these en- 
tropies with appropriate constraints, i.e., the normalization of distribution function: 
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J daP^^"^ (a) = const, and the g-variance being kept constant as a known quantity: 
a2 = {{a - aof)q = (/ daP^i) - aof)/^daP(^\aY, we have 

pW(a)oc[l-(a-ao)7(Zi«)2]^/(^-') (2) 

with {AaY = 2X/[(1 — (jf)ln2]. The Renyi entropy has an information theoretical 
basis, and has the extensive character as the usual thermodynamical entropy does. 
On the other hand, the Tsallis entropy is non-extensive, and therefore it provides 
us with one of the attractive trials for generalized statistical mechanics which deals 
with, for example, systems with long range correlations with a hierarchical structure 
where the usual extensive characteristics in statistical mechanics are not applicable. 
Note that the distribution functions which give an extremum of each entropy have a 
common structure (see (H) for the present system where q < 1) in spite of different 
characteristics of these entropies, and that the values of a are restricted within the 
range [amin, ctmax], where a^ax - ao = ao - amin = Aa. Note that = 2X/[(3 - 
g)ln2]. 

By making use of an observed value of the intermittency exponent /i as an input, 
the quantities ao, X and the index q can be determined, self-consistently, with the 
help of the three independent equations, i.e., the energy conservation: (e„) — e, the 
definition of the intermittency exponent ji: (e^) = e^5~^ , and the scaling relation: 
1/(1 — 9) = 1/q;_ — l/a+ with a± satisfying f{a±) = p^-p^. Here, the average (• • •) 
is taken with P^") {a) . The scaling relation is a generalization of the one derived first in 
p^ , p9[ | to the case where the multifractal spectrum has negative values. For the region 
where the value of fj, is usually observed, i.e., 0.13 < /i < 0.40, the three self-consistent 
equations are solved to give the approximate equations ao — 0.9989 -I- 0.5814/i, X — 
-2.848X 10-3+1.198/1 and g = -1.507+20.58/i-97.11/i^-l-260.4/i3_365.4Ai4-|-208.3Ai5. 
These equations are slightly different from those given in |p^ , since the region of fi 
has been extended a little bit. 




FIG. 1. The longitudinal scaling exponents (^l^ of the velocity structure function. The 
second figure (b) is a magnification of the first (a). The results of DNS obtained by Gotoh 
et al. {R\ — 381) are shown by closed circles. The present theoretical formula was used 
to determine the value fi = 0.240 of the intermittency exponent, and the result is drawn by 
solid line. Dotted line represents K41, whereas dashed line She-Leveque. The prediction of 
the log-normal model is given by short-dashed line with /i = 0.240. 



If we put q = 1 from the beginning, i.e., starting with the Boltzmann-Shannon 
entropy, and take an extremum with the same constraints given above for q = 1, we 
have a Gaussian distribution function. The parameters ao and X are then determined 
by two conditions, i.e., energy conservation and the definition of the intermittency 
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exponent. The obtained results completely coincide with those in the log- normal 
model @-|^. Therefore, the present approach can be interpreted as an extension of 
the one in the log-normal model. 
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FIG. 2. The transverse scaling exponents ("^ of the velocity structure function. The second 
figure (b) is a magnification of the first (a). The results of DNS obtained by Gotoh et al. 
(-Ra ~ 381) are shown by closed circles. The present theoretical formula (a) was used to 
determine the value ^ — 0.327 of the intermittency exponent, and the result is drawn by 
solid line. Dotted line represents K41, whereas dashed line She-Leveque. The prediction of 
the log-normal model is given by short-dashed line with /i = 0.327. 



It is known that in turbulent flow there are two mechanisms to rule its dissipative 
evolution, i.e., the one controlled by the kinematic viscosity that takes care thermal 
fluctuations, and the other by the turbulent viscosity that is responsible for intermit- 
tent fluctuations related to the singularities in velocity gradient. Therefore, it may 
be reasonable to assume that the probability 77'^"^(x„)(ix„ to find the scaled velocity 
fluctuation — Su„/Suq in the range x„ ^ a;„ + dxn can be divided into two parts: 

n^^\xn)dXn = n^\xn)dXn + 77s'"\ |.T„ | . (3) 

Here, the normal part PDF n^^\xn) stemmed from thermal dissipation, and the 

singular part PDF 77g"-'(|a;„|) from multifractal distribution of the singularities. The 

latter is derived through 77g"''(|a;„|)(ix„ = P''-"-'^ {a)da with the transformation of the 

variables: \xn\ = S"^^ ■ 

The mth moments of the velocity fluctuations, defined by ((|a;n|"')) — 
/^^2;„|a;„|'"77(")(a;„), are given by 

((KI")>-27(:) + (l-27<"V™e (4) 
where 03^ = {2/[Cl^^{l + 0^^)]}^^^ with Cq = 1 + 2f{l - q)X\n2, and 

/OC 
dXn \Xnrn^"\x,,). (5) 
-00 

We used the normalization: ((1)) = 1. The quantity 
— Qoto/S — 2Xrri^ / 

is the so-called scaling exponent of the velocity structure function, whose expression 
was derived first by the present authors ||l5|-|l8||. In this paper, we will use the formula 



9(l+^S - l-log2(l + Ci/;3)]/(l-g) (6) 
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in order to extract the value of the intermittency exponent fi for the best fit to 
the measured scaling exponents by the method of least squares. 




FIG. 3. Experimentally measured PDF of the longitudinal velocity fluctuations by Gotoh 
et al. for Rx = 381 are compared with the present theoretical results fl^"^ Closed circles 
are the symmetrized points obtained by taking averages of the left and the right hand sides 
DNS data. For the experimental data, the separations r/rj — in/rj are, from top to bottom: 
2.38, 4.76, 9.52, 19.0, 38.1, 76.2, 152, 305, 609 and 1220. Solid lines represent the curves 
given by the present theory with q = 0.391 (/i = 0.240). For the theoretical curves, the 
number of steps in the cascade n are, from top to bottom: 21.5, 20.0, 16.8, 14.0, 11.8, 10.1, 
9.30, 8.10, 7.00 and 6.00. For better visibility, each PDF is shifted by -1 unit along the 
vertical axis. 



There usually appears an asymmetry (a negative skewness) in the PDF's of the 
longitudinal velocity fluctuations, but this is not so vivid in the PDF's of the transverse 
ones (see, for example. Fig. 15 and Fig. 16 in |^). Assuming that the singularities 
in the velocity gradient may contribute mainly to the symmetric part, we will deal 
with the symmetric part of the PDF's of velocity fluctuations in the following. 

With the help of the new variable 

^„ = suj{{6uir/' = xjiixir/' = (7) 

scaled by the variance of velocity fluctuations, the PDF -/7^"-'(|^n|), introduced through 



^^"^(l^n|)d^„ = ir(")(|x„|)da;„, is given by p 20| 



fi^'^H^n) ^ n^lHu forie„i<e* (s) 

n'^^Hin) = nl<\Cn) for C < \Cn\ < Un-'^'^''\ (9) 
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Here, ^„ — [272"'' (5^''^ + (1 — 27Q"^)a2]~"'^/^. Assuming that, for smaller velocity 
fluctuations |^„| % the contribution to the PDF of the velocity fluctuations comes 
mainly from thermal fluctuations related to the kinematic viscosity, we take for the 
PDF n'^^}{in) a Gaussian function i.e., 

= 77s(")e-[W("*)l[(«"/«;)'-il/2 (10) 

with Tig'"' = 3(1 - 27^"V(2^n\/27r^|ln(5n|). On the other hand, we assume that 
the main contribution to tl^iXn) may come from the multifractal distribution of 

n/(l-q) 



singularities [p^-^ related to the turbulent viscosity, i.e., fli^i^n) — ^s"^(|6i|)- 



i^)fC \ — fT(") ^" 



31n|gn/6K0 
Aa I In Sn I 



(11) 



with |^„,o| = CnSn°^^ The point ^* was defined by ^* = ^„(5n ^'^ ^^^^ where 

a* is the solution of C2/2 - a/3 + 1 - /(a) = 0, and i7^"'(C„) and /li<'(C„) were 
connected at ^* under the condition that they should have the same value and the 
same derivative there. Here, 

/(a) = 1 + (1 - q)'' log, [l-{a~ a^f / {Aaf] (12) 

is the multifractal spectrum derived by the relation p(")(a) oc Sn'^'-"^ §,0, 

that reveals how densely each singularity, labeled by a, fills physical space. 

Now, we proceed to analyze the data of DNS conducted by Gotoh et al. |^l| for fully 
developed turbulence. In the following theoretical analyses, we will take 6 = 2 for the 
number of pieces of "eddies" generated at each step in the energy cascade. In Fig. |l| 
and Fig. ^ we put, respectively, the measured scaling exponents for longitudinal 
velocity fluctuations and Cm for transverse ones at Ra = 381 (closed circles) Q. The 
scaling exponents (^) derived by the present theory are given by a solid line in each 
flgure. There are also represented, as references, the predictions of K41 (dotted line) 
1], of She-Leveque (dashed line) and of the log- normal model (short-dashed line) 
4--|^. We determine, self-consistently, the values of the intermittency exponent fj, by 
fltting the formula (^ with the ten observed values of the scaling exponents, Cm and 
Cm I with the help of the method of least squares. The determined values of /i are 
listed in Table ^ with the values of the corresponding parameters q, and X . The 
latter three parameters are obtained by the /i-dependent functions that are given in 
this paper as the solutions of the self-consistent equations. Note that the theoretical 
lines thus determined locate within the experimental error bars at each ten observed 
points. Therefore, the relation /i = 2 — Ce is satisfied within the experimental error 
bars. We obtain a+ - ao = ao - a- = 0.6818 (0.8167), Aa = 1.160 (1.566) for the 
longitudinal (transverse) fluctuations. 

There is an argument that the scaling exponents C2 l^r the longitudinal velocity 
structure function and Cj for the transverse function should be equal for the isotropic 
and incompressible turloulence. In Table |, we see that C2 = 0.696 and Cj = 0.707 
within the present theoretical analysis, which give, respectively, 1.696 and 1.707 for 
the exponent of the Kolmogorov spectrum which is 5/3 = 1.6 in K41 ||l|. The small 
deviation (1%) between C2 and Cj can be attributed to the finite sample size and 
the small amount of fiow anisotropy [^T| . The scaling exponent C4 of the fourth order 
moment, related to the second order pressure structure function, is reported in [ pl| as 
C4 = 1.30±0.02 and Cj 1.26±0.03 at Rx = 381. These values are comparable with 
the ones in Table |, i.e., C4 = 1.277 and Cj = 1.256 derived by the present theory. 
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FIG. 4. Experimentally measured PDF of the transverse velocity fluctuations by Gotoh et 
al. for Rx = 381 are compared with the present theoretical results n^"\(^„). Closed circles 
are the symmetrized points. The separations r/r; = in/v are the same as in Fig. |^. Solid 
lines represent the curves given by the present theory with q — 0.543 {fi — 0.327). The 
number of steps in the cascade n are, from top to bottom: 20.0, 18.2, 14.6, 11.4, 9.30, 7.90, 
6.80, 5.60, 4.70 and 4.00. 

The comparison between measured PDF's of the velocity fluctuations at R\ = 381 
by DNS ||2l[ and those obtained by the present analysis jTsj-Q is given in Fig. ^ 
for longitudinal fluctuations and in Fig. |^ for transverse ones. In order to extract 
the symmetrical part from the DNS data, we took mean average of those on the 
left hand side and on the right hand side. The symmetrized data are described by 
closed circles. The solid lines are the curves n^"\^n) given by (|lO| ) and ( |ll| ) with the 
values of parameters in Table ^ The values of ^* are about 1 to 1.5. The separations 
r/r] = Inh are, from top to bottom: 2.38, 4.76, 9.52, 19.0, 38.1, 76.2, 152, 305, 609 
and 1220 for DNS |2^ data. On the other hand, the number n of steps in the cascade 
for longitudinal (transverse) fluctuations are, from top to bottom: 21.5, 20.0, 16.8, 
14.0, 11.8, 10.1, 9.30, 8.10, 7.00 and 6.00 (20.0, 18.2, 14.6, 11.4, 9.30, 7.90, 6.80, 5.60, 
4.70 and 4.00). These values are obtained by the method of least squares with respect 
to the logarithm of PDF's for the best fit of our theoretical formulae ( |l0| ) and ([Tl| ) 
to the observed values of the PDF's by discarding those points which have observed 
values less than 10^^ (10~*) since they scatter largely in the logarithmic scale. We 
see an excellent agreement between the measured PDF's and the analytical formula 
of PDF derived by the present self-consistent theory. 

The dependence of n on r/r] for longitudinal (transverse) fluctuations, extracted 
from Fig. ^ (Fig. U), is shown in Fig. ^ (a) (Fig. |^ (b)) by solid and dashed lines. 
These lines are obtained by the method of least squares within linear fits, and are 
given by 
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n = -1.050 X log2 r/rj + 16.74 (for < r), 
n = -2.540 X log2 r/r] + 25.08 (for r < i^) 



(13) 
(14) 



with the crossover length £^/i] = 48.26 for longitudinal fluctuations, and 

n = -0.9896 xlog2r/?7 + 13.95 (for < r), (15) 
n = -2.820 X log2 r/77 + 23.87 (for r < ej) (16) 

with the crossover length £]: /rj — 42.57 for transverse fluctuations. We see that 
/r] ^ i];/r], and that these crossover lengths have values close to the reported value 
of the Taylor microscale A/77 = 38.33 in at R\ = 381. The value of ij /i] is 
very close to the one in Fig. 31 of (21]. Whereas, the values of i^/rj is different 
from the one obtained in ||2l| where Gotoh et al. got the relation i^/rj w 2£]:/ri 
within their analyses. In spite of the difference in the values of t]^/rj^ the inequality 
/r] < £]:/ rj may be attributed to a manifestation of the incompressible nature of the 
fluid under consideration as interpreted in . Notice that this kind of crossover was 
not observed in the analysis of the experiment conducted by Lewis and Swinney 
for turbulent Couette- Taylor flow. 




10" 



10" 



r/r] 




FIG. 5. Dependence of n on r/r; is given by closed circles, 
for longitudinal fluctuations. The crossover occurs at (^/rj 
for transverse fluctuations. The crossover occurs at /ri = 42.57. 



r/r) 

(a) Those extracted from Fig. 3 
= 48.26. (b) Those from Fig. 4 



The equations (|l3| ) and ( |l5| ) have the slope close to —1 with respect to log2r/?7 
that is consistent with (^ for 8 — 2. Therefore, we conclude that the upper scaling 
region £]^/r] < r/rj % 1220 /rj < r/rj < 1220) found within the present systematic 
analysis can be interpreted as the inertial range for longitudinal (transverse) fluctu- 
ations. Note that this inertial range for longitudinal fluctuations is wider than the 
longitudinal scaling range 2A/?7 < r/r] < 220.7 observed in |2^ at R\ = 381. The 
latter estimate was done by looking for the flat region of the observed Cm with re- 
spect to the separation r/rj. The existence of the crossover length was observed first 
by Gotoh et al. in the analyses of their DNS data [|l| as the low end of the scaling 
range for the longitudinal or the transverse fluctuations, and was attributed to the 
existence of structure of this order in turbulence, e.g., a shear layer or a vortex tube. 

It is remarkable that in Fig. ^ we can see, within the present self-consistent analysis, 
the existence of another new scaling region, the lower scaling region, i.e., 2 r/rj < 
i^/f] (2 r/rj < ij /rj) for longitudinal (transverse) velocity fluctuations. It may be 
interesting to point out that the slope in (|l^) and (|l^) can be —1 with respect to 
log^r/rj with S — 5^ = 1.33 for longitudinal fluctuations and with S = S"^ = 1.28 
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for transverse ones, respectively. Therefore, in this lower scahng region, it can be 
interpreted that eddies break up, effectively, into 6^ {6^) pieces at each step of the 
energy cascade for longitudinal (transverse) fluctuations. 

Assuming that (|^) with (5 = 2 is applicable in the "inertial range" (^/ri,£]:/ri < 
r/r] < io/v^ we can extract from (|l^) and (|l^) the values of the Reynolds number 
Re^ = 5.236 x 10^ for longitudinal eddies, and of Re""" = 7.113 x 10^ for transverse 
eddies. If we adopt the relation Rl = A^Re^ = A^Re^ Hi, we have A^ = 0.02772 
and A'^ = 0.2041. 

In this paper, we showed how precisely the formulae, derived by our theory [p"5|-p0|, 
can explain the data observed in the DNS |^l[]. We found, explicitly, that there exist 
two scaling regions, i.e., the upper scaling region with larger separations which may 
correspond to the scaling range observed by Gotoh et al. , and the lower scaling region 
with smaller separations which is a new scaling region extracted first by the present 
systematic analysis. These scaling regions are divided by a definite length approxi- 
mately of the order of the Taylor microscale A, that may correspond to the crossover 
length introduced by Gotoh et al. |21 1 as the length at which a scaling behavior of the 
structure functions ceases due to the effect of kinematic viscosity. We also saw that 
the inertial range (the upper scaling region) for longitudinal fluctuations is wider than 
the scaling range extracted by Gotoh et al. The discovery of the lower scaling 

region by the present analysis may tell us that the distribution of singularities of the 
velocity gradient in physical space is robust enough to produce a scaling behavior 
even in the dissipation range. In other words, the system of turbulence may have an 
intrinsic scaling property, applicable in a wider region than the inertial range, whose 
origin can be attributed to the multifractal distribution of the singularities in the ve- 
locity gradient. A study of the PDF's for velocity derivatives by the present analysis 
is one of the attractive open problems and is now in progress. It is also one of the 
attractive future problems to provide with a dynamical approach based on the aspect 
given by the present ensemble theoretical approach in order to interpret, especially, 
the lower scaling region. These will be reported elsewhere. 

The authors are grateful to Prof. T. Gotoh for enlightening discussion and his 
kindness to show his data prior to publication. The authors would like to thank 
Prof. C. Tsallis and Dr. A.K. Rajagopal for their fruitful comments with encourage- 
ment. 
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longitudinal 


transverse 


A* 


0.240 




0.327 




<? 


0.391 




0.543 






1.138 




1.189 




X 


0.285 




0.388 




m 


Sm 




Sm 






present theory 


DNS 


present theory 


DNS 


1 


0.3637 


0.370±0.004 


0.3747 


0.369±0.004 


2 


0.6965 


0.709±0.009 


0.7073 


0.701±0.01 


3 


1.000 


1.02±0.02 


1.000 


0.998±0.02 


4 


1.277 


1.30±0.02 


1.256 


1.26±0.03 


5 


1.529 


1.56±0.03 


1.480 


1.49±0.04 


6 


1.761 


1.79±0.04 


1.674 


1.69±0.05 


7 


1.973 


1.99±0.04 


1.843 


1.86±0.05 


8 


2.169 


2.18±0.04 


1.990 


2.00±0.04 


9 


2.350 


2.35±0.04 


2.117 


2.11±0.05 


10 


2.519 


2.49±0.04 


2.227 


2.20±0.06 



TABLE L Comparison of the values of the scaling exponents (|6|) derived by the present 
theory and those for the longitudinal and transverse velocity fluctuations observed in DNS 
conducted by Gotoh et al.. The value p for each fluctuations is also listed with the corre- 
sponding values for g, ao and X. 
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